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A set of general constructing schemes are unveiled to predict a large family of stable boron monoelemental, 
hollow fullerenes with magic numbers 32 + 8/c {k > 0). The remarkable stabilities of these new boron fullerenes 
are then studied by intense ab initio calculations. An electron counting rule as well as an isolated hollow rule 
are proposed to readily show the high stability and the electronic bonding property, which are also revealed 
applicable to a number of newly predicted boron sheets and nanotubes. 
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Carbon nanostmctures such as fullerenes, nanotubes and 
graphenes Mia have been extensively studied in the 
past decades. As boron and carbon are neighbors in the peri- 
odic table, and possess many structural analogies i^lzl], it is 
expected that boron could also form nanostructures similar to 
carbon. The boron clusters, nanotubes (BNT) and sheets are 
thus actively explored in recent years ||8l Isi llOfl. More cur- 
rently, a stable boron fullerene Bgo was predicted 111 ill, and 
a new type of boron sheet (NBS) and related BNTs |tl2l ll3|] 



were found to be remarkably more stable than the boron sheet 
and BNTs with triangular structures; several other fullerene- 
like boron nanostructures are also predicted In this 

Letter, we shall present generic constructing schemes that can 
produce a large family of novel stable boron nanostructures 
that include Bgo buckyball and the NBS as members. 

Let us start with analyzing the geometrical structures of Bgo 
and the NBS. They contain three basic motifs: hollow pen- 
tagon (HP), hollow hexagon (HH), and filled hexagon (FH) 
(i.e. hexagon with an additional atom in center), respectively. 
A closer inspection reveals that both Bgo and NBS are entirely 
composed of such a snowdrop-like structure [Fig. 1(a)] in that 
a central FH is surrounded by three FHs and three "hollows" 
(HPs or HHs), and every FH can be viewed as the center of 
the structure. For Bgo and the NBS, the HPs and HHs act as 
"hollows", respectively [Figs. 1(b) and (c)]. In light of this 
observation, we shall show that such a snowdrop-like struc- 
ture plays a critical role in the stability of Bgo and NBS as 
well as other boron quasi-planar structures (fullerenes, sheets, 
nanotubes, etc), and thus can be utilized as building blocks to 
construct novel boron nanostructures. Henceforth, the boron 
fullerenes that are composed of snowdrop-like structures are 
coined as S -fullerenes for brevity. 

Geometrically, in an S-fullerene, every HP should be sur- 
rounded by five FHs, and every HH is surrounded by six FHs, 
whereas there are three "hollows" (HPs or HHs) around every 
FH. Thus we have 5Nhp + QNhh = ^NpH, where Nhp, 
Nhh, and NpH denote the number of HP, HH, and FH, re- 
spectively. If we take the S-fullerene as a polyhedron, and 
each HP, HH or FH as a surface, then Euler theorem, e (edges) 
+ 2 = V (vertices) + / (faces), leads to Nhp = 12, implying 
that in an S-fullerene the number of hollow pentagons should 
be exactly 12. Consequently, we have NpH = 20 + 2Nhh- 




FIG. 1: (Color online) The snowdrop-like structures in Bso and 
new boron sheet (NBS). (a) a snowdrop-like motif, in which a filled 
hexagon (FH) is surrounded by three FHs and three "hollows" (pen- 
tagons or hexagons, denoted by circles); (b) Bso; (c) NBS, the yellow 
parallelogram indicates an unit cell of the NBS. The blue atoms form 
the snowdrop-like units in Bso and NBS. 




FIG 2: (Color online) Illustration of leapfrog and modified leapfrog 
transformations. (a) a three-connected vertex becomes a new 
hexagon (in leapfrog) or a hexagon with an additional atom in cen- 
ter (in modified leapfrog); (b) a C20 converts to Ceo and Bso by the 
leapfrog and modified leapfrog transformations, respectively; (c) a 
C24 transforms to C72 and Bge. 



In an S-fullerene the total number of boron atoms should be 
n = [bNHP + Q{Nhh + Nfh)]/S + Nfh = 80 + SNhh- 
For B80, since Nhh = 0, Nfh = 20, and n = 80. 

These relations put restrictions in searching for S- 
fullerenes, but they are still inadequate to find the concrete ge- 
ometrical structures. Fortunately, through intensive attempts 
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we disclose that a modified leapfrog transformation (MLT) is 
useful to generate S-fullerenes. The original leapfrog transfor- 
mation lll6l] has been applied to construct the so-called elec- 
tronic closed-shell fullerenes Ceo+e/c (leapfrog-fullerenes). 
To produce S-fullerenes, the only thing we should do is to 
modify the operations shown in Fig. 2(a). In the original 
leapfrog operations, every vertex transforms to an HH (con- 
taining six atoms). Here we propose an MLT, in which every 
vertex converts to an FH (containing seven atoms), suggesting 
that every vertex of the original polyhedron provides an extra 
atom. Since these extra atoms could compose a polyhedron 
identical to the original one, the MLT can be considered as the 
combination of the original polyhedron and its leapfrog. For 
instance, Bgo, the modified leapfrog of C20, could be viewed 
as a combination of C20 and its leapfrog Ceo structures [Fig. 
2(b)] . As a result, in the final modified leapfrog structure, the 
number of atoms grows four times vertices as many as the 
original, i.e., 80 + 8/c (/c > 0, /c 7^ 1), while the symmetry 
is still kept. The most important point is that all structures 
generated by the MLT fall into the S-fullerenes. Through a 
reverse transformation, every S-fullerene could be converted 
to the original polyhedron. Therefore, Bgo+s/c (k > O^k 1) 
can be produced from the structure of C2o+2fe (k > O^k 1) 
by means of the MLT [Figs. 2 (b) and (c)]. The isomer num- 
bers of both Ceo+e/c and Bgo+s/e are equal to that of C20+2/C 
for each k (k > 0,/c 7^ 1), and the numbers of valence elec- 
trons are identical for Ceo+efe and Bgo+s/c with the same k. 
In addition, the MLT may also be applied to generate other 
quasi-planar structures composed of snowdrop-like units. As 
the leapfrog of a graphene structure is still graphene, the MLT 
of graphene leads to the NBS structure; similarly, the struc- 
tures of BNTs can be generated by the MLT of CNTs. 

Several structures of boron S-fullerenes as examples shown 
in Fig. 3 have been generated and then studied by means of ab 
initio calculations [17]. The isomer numbers and symmetries 
of Bgo+sfe are listed in Table I. For /c = 6, 7, 8, 9, 10, the 
isomer number is 6, 6, 15, 17, 40..., respectively, which grows 
rapidly with increasing k. For instance, when k = 0, the cor- 
responding S-fullerene is the known spherical Bgo 111111 : when 
k = 2, it produces Bge [Fig. 3(a)], which bears Ded symmetry 
and has a round pillow shape with two HHs located at the De^ 
axis; when /c = 3, it gives a pineapple shaped B104 with D^h 
symmetry [Fig. 3(b)]; when /c = 4, it generates B112 [Figs. 
3(c)-(d)], which has two isomers: one possesses D2 symmetry 
with an elliptical pillow shape, while another has T^^ symme- 
try and looks like a regular tetrahedron with four smoothed 
vertices; when k = 5, it yields B120 [Figs. 3(e)-(g)], which 
has three isomers with one bearing D^h symmetry, and the 
other two having C2v symmetry. The D^h B120 has a cylin- 
der capsulate shape, and can be viewed as a short open-end 
BNT B40 unit (which is rolled from NBS) capped with two 
halves of Bgo- This capped BNT can be elongated by repeat- 
ing the middle open-end BNT B40 unit to generate similar 
D^h symmetry capped BNTs B8o+40t (t > 1) CI [3, all of 
which belong to a subset of Bgo+sfe- Generally, using a proper 
S-fullerene as the cap, and the open-end BNT with different 
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FIG. 3: (Color online) Schematic structures and Mulliken atomic 
charge populations of several boron S-fullerenes. The unit of color 
bar is electrons. 



TABLE I: The number (N^so), symmetries, relative energies per atom 
(E, in eV, relative to that of graphene and NBS, respectively), and 
HOMO-LUMO gap energies ( AE, in eV) of Ceo+efc and S-fullerenes 
Bso+sfc (k<5) isomers (obtained with Gaussian03). 
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20 



5 3 



D2 

D5/^ 



72 0.365 
78 0.346 
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0.141 0.927 

0.137 0.709 

0.136 0.828 

0.135 0.555 
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sizes as the body, the other similar subsets of boron nanotubes 
could be generated. 

The energies of Ceo+e/c and Bgo+sfc (k 1, k < b) rela- 
tive to graphene and NBS, respectively, are listed in Table I. 
Note that different isomers with the same k may have differ- 
ent structures and symmetries. It can be seen that the rela- 
tive energies decrease with the increase of k for both Ceo+efe 
and B80+8/C, which are expected to converge to the values of 
graphene [16] and NBS, respectively, suggesting that the sta- 
bilities of both fullerenes increase with the growth of size. 
Thus, Bgo+s/c (k > 0, k 1) gives rise to a large class of 
stable boron fullerenes, where Bgo is just the first member 
of this family, while the NBS is the infinite k analogue. To 
check further the stabilities of these boron fullerenes, we cal- 
culate the total energies per atom for the caged Bgo+s/c (k y^ I, 
k < b) with other different structures that do not belong to the 
S-fullerene family, and find that they are all less stable than 
those of S-fullerenes. Note that Bge is about 0.019 eV/atom 
more stable than Bgo in energy, and the members with larger 
k have even more lower relative energies. The present class of 
S-fullerenes are more stable than B180 proposed earlier [14]. 

In the preceding S-fullerenes, only are HP and HH con- 
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sidered as boron "hollows". If we take the hollow quad- 
rangle (HQ) into account, then the S-fullerenes would be 
extended. For the boron fullerenes composed of HQ, HP, 
HH, and FH, the total number of boron atoms satisfies n = 
[4^Nhq + ^Nhp + 6{Nhh + Nfh)]/3 + Nfh, while Euler 
theorem requires Nhp + "^Nhq — 12, and Npu = 20 + 
'1{Nhh - Nhq). Then, we have n = 80 + S{Nhh - Nhq). 
Consider a few special cases, (a) If Nhq = 0, it recovers to 
the preceding Bgo+s/e (k > 0, k 1); (b) If Nhp = 0, we 
get Nhq = 6, Nfh = 8 + 2Nhh, and n = 32 + 8Nhh, 
which leads to 632+8/0 (k > 0) that consist of HQs, HHs 
and FHs; (c) If Nhh = 0, we have Nfh = 20 - 2Nhq, 
n = 80 — 8Nhq, and < Nhq < 6, which corresponds to 
B32+8/C (0 < A: < 6) that are composed of HQs, HPs and FHs. 
As examples. Oh B32, B40, D2d B48, D3/, B56, B64, 
and Dsh B72 can be obtained, as shown in Figs. 3(i)-(l); (d) If 
Nhh — Nhq = 1, then B88 can be formed, which has abun- 
dant isomers because Nhh and Nhq can have various com- 
binations to satisfy such a difference of unity. The structure 
of one C2v B88 isomer is given in Fig. 3(h). Ab initio calcula- 
tions show that these several isomers also follow the trend that 
relative energies per boron atom decrease with the increase of 
k. Thus if quadrangles were added to the set of "hollows", we 
would have a class of fullerenes B32+8fc (k > 0), and their 
isomers also increase largely for every k, which can be gen- 
erated by the MLT as well, and enriches the family of boron 
fullerenes. 

By carefully checking the structures of the S-fullerenes and 
NBS, we uncover an electron counting rule (ECR) and an iso- 
lated hollow rule (IHR) that would manifest the fundamen- 
tal importance of the snowdrop-like structure and, reasonably 
explain the high stability and reveal the electronic bonding 
property of these boron nanostructures. Let us look again at 
the structures of B80 and NBS. The former is composed of 12 
HPs and 20 FHs, while a unit cell of NBS has one HH and 
two FH [Fig. 1(c)]. As one FH can be viewed as a group 
of six triangles, B80 contains 120 triangles and an NBS unit 
cell has 12 triangles. Since each boron atom has three va- 
lence electrons, one may find that the total number of valence 
electrons of B80 is 240, just twice the number of triangles. In- 
terestingly, the NBS has the same electron counting feature, 
as one unit cell has 8 boron atoms, contributing 24 valence 
electrons, which is just also twice the number of triangles it 
possesses. It is easy to substantiate that all of the S-fullerenes 
follow this electron counting scheme. In the central FH of 
a boron snowdrop-like unit, the atoms on six periphery ver- 
tices are shared by the neighboring FHs, and only is the cen- 
tral atom monopolized, so the net number of boron atoms that 
the central FH has is only four, and the corresponding num- 
ber of valence electrons is 12, just twice the six triangles it 
possesses. Since all FHs in the S-fullerenes can be viewed 
as the central FHs of snowdrop-like structures, being equiva- 
lent to the above electron counting scheme, the total number 
of valence electrons of S-fullerenes are twice the number of 
triangles they have. 

Such an agreement on electron counting property is not ac- 



cidental, which could be rationalized in a simple way 
As is well-known, the multi-center deficient-electron bond- 
ing scheme exists widely in boron structures ||2Q|] . Suppose 
that three boron atoms on the vertices of a triangle share a 
three-center two-electron bond, i.e., a triangle would consume 
two electrons, then the above ECR could be properly under- 
stood in a way that, in stable quasi-planar boron structures 
that are composed of FHs, HPs and HHs, etc., the total num- 
ber of the valence electrons is equal to twice the number of 
triangles they contain, leading to that the number of boron 
atoms is as quadruple as the number of FHs (n = 4:Nfh)- 
The ECR is thus consistent with the snowdrop-like structure, 
as the latter is a kind of the geometrical realization of ECR, 
whilst the ECR can be considered as the electronic interpre- 
tation of the snowdrop-like unit. On the other hand, in all 
S-fullerenes, the "hollows" (HPs and HHs) are separated by 
FHs, as the requirement that every FH should be the center of 
a snowdrop-like structure cannot be satisfied if there exist any 
adjacent "hollows". Therefore, an IHR for boron S-fullerenes 
can be expected. The leapfrog procedure never generates ad- 
jacent pentagons [16], and similarly, the MLT never produces 
adjacent "hollows". It is interesting to note that the NBS 
and the related BNTs also follow the ECR and IHR, which 
suggests that the above two rules may be applicable to other 
boron nanostructures. While the IHR calls for boron "hollow" 
be well- separated by FHs, the ECR requires a quantitative 
constraint between the number of boron "hollows" and that 
of FHs. Consequently, the IHR and ECR are mutually con- 
strained, which results in appropriate boron nanostructures. 

The electronic structures of B8o+8fc are obtained by the ab 
initio calculations. Figure 3 shows the Mulliken atomic charge 
populations of several S-fullerenes. The boron atoms in cen- 
ters of FHs are marked in green, showing they are positively 
charged, while most of others are red, indicating they are neg- 
atively charged. Recall that all carbon atoms in Ceo+efe are 
electroneutral. With regard to electron transfers in B8o+8/c, 
the central atoms of FHs act as electron donors, and the other 
atoms act as electron acceptors, which are consistent with the 
donor- acceptor hypothesis on NBS 1IT2I1 . and also the charge 
transfer in monoatomic systems 1I21I1 . 

The gap energies between the lowest-unoccupied- 
molecular-orbital (LUMO) and the highest-occupied- 
molecular-orbital (HOMO) of Ceo+e/c and B8o+8/c 
(0 < k < 5, k 7^ 1) are listed in Table I. Interestingly, 
the isomers of both fullerenes with k = 0(/^), 4(Tc^), ^{C^) 
have larger gap energies than their near neighbors. Fig. 4(a) 
shows the energy levels around the HOMO-LUMO gaps of 
Ceo+e/c and B8o+8/c (k = 0, 2, 3). For the major energy levels 
of B8o+8fc, one may find correspondences to their Ceo+e/c 
counterparts with the same degeneracy and symmetry. For 
instance. Ceo and B80 both have a five-fold degenerate hu 
HOMO and a three-fold degenerate hu LUMO [22]; C72 
and B96, C78 and B104 also exhibit similar corresponding 
relations as indicated in Fig. 4(a). The visualizations of the 
HOMO and LUMO of B8o+8fe (k = 0, 2, 3) are presented in 
Figs. 4(b)-(d). We find the main profiles of corresponding 
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FIG. 4: (Color online) (a) Energy levels near HOMO-LUMO gaps 
of Ceo, Bso, C72, B96, C78, and B104; The profiles of HOMO and 
LUMO orbitals of (b) Bso; (c) Bge and (d) B104. 



MOs (Ceo+e/c MOs are not shown here) exhibit striking 
similarity between Ceo and Ego, C72 and Bge, C78 and B104, 
and so forth. We believe that such an observation from the 
first several members of Bgo+s/c family would be kept for the 
members with even larger k. Since the frontier MOs dom- 
inate the primary chemical property of a molecule, Bgo+s/c 
may have chemical properties similar to the corresponding 
Ceo+e/c- Besides, as the carbon fullerenes can be stabilized 
by the exohedral chemical deriving and metallic endohedral 
method i3 Q Us], the stability of boron fullerenes might 
also be improved in the same manner. 

In summary, the general constructing schemes were pro- 
posed to construct a large family of stable boron fullerenes, 
whose remarkable stabilities were confirmed by intensive 
first-principles simulations. The empirical ECR and IHR were 
suggested to readily explain the geometrical stability and re- 
veal the electronic bonding of the predicted boron fullerenes, 
which can also be shown applicable to boron quasi-planar 
structures such as nanotubes and sheets, and be useful for 
seeking for other stable boron nanostructures. The present 
predictions are awaited to synthesize experimentally, which 
may bring about a new dimension of boron chemistry. 

All of the calculations are completed on the supercomputer 
NOVASCALE 6800 in Virtual Laboratory of Computational 
Chemistry, Supercomputing Center of Chinese Academy of 
Sciences. This work is supported in part by the National Sci- 
ence Foundation of China (Grant No. 10625419). 
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